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ABSTRACT 

We employ numerical simulations and simple analytical estimates to argue that dark matter substructures 
orbiting in the inner regions of the Galaxy can be efficiently destroyed by disk shocking, a dynamical process 
known to affect globular star clusters. We carry out a set of fiducial high-resolution collisionless simulations 
in which we adiabatically grow a disk, allowing us to examine the impact of the disk on the substructure 
abundance. We also track the orbits of dark matter satellites in the high-resolution Aquarius simulations and 
analytically estimate the cumulative halo and disk shocking effect. Our calculations indicate that the presence 
of a disk with only 10% of the total Milky Way mass can significantly alter the mass function of substructures in 
the inner parts of halos. This has important implications especially for the relatively small number of satellites 
seen within ^30 kpc of the Milky Way center, where disk shocking is expected to reduce the substructure 
abundance by a factor of 2 at 10 9 M and 3 at 10 7 M Q . The most massive subhalos with 10 10 M Q survive 
even in the presence of the disk. This suggests that there is no inner missing satellite problem, and calls into 
question whether these substructures can produce transient features in disks, like multi-armed spiral patterns. 
Also, the depletion of dark matter substructures through shocking on the baryonic structures of the disk and 
central bulge may aggravate the problem to fully account for the observed flux anomalies in gravitational lens 
systems, and significantly reduces the dark matter annihilation signal expected from nearby substructures in 
the inner halo. 

Subject headings: Galaxy: disk - Galaxy: formation - Galaxy - dynamics - Galaxy: structure 



1. INTRODUCTION 

In the cold dark matter (CDM) scenario, structure grows 
hierarchically, with small objects collapsing fir st and continu- 
ously merging to form larger and larger bodies (IWhite & Reesl 
1978). Numerical simulations have tracked the evolution of 
this process, beginning at early times with objects havin g 
masses comparable to that of the Earth ( Diemand et al.l2005h . 
and progressing through many orders of magnitude in mass 
to the scales of gal axies, galaxy clusters, and cosmic large - 
scale structure (e.g. iGhigna et alJl2000t ISpringel et al.ll2005l) . 
In the Milky Way, the cumulative number TV of dark matter 
satellites of a given mass M sat is predicted to scale with the 
subhalo mass as N<xM~°- 9 (Sprin gel et alll2008bl) . Because 
the subhalos hosting the satellite galaxies formed early, when 
the Universe was dense, the smaller structures are thought 
to be resilient to tidal disruption and the simulations pre- 
dict that ^100 with maximum circular velocity larger than 
20kms~' should surv ive to the present day for a halo like 
that of the Milky Way dKlypin et aljl 19991: [Moore et alj|1999t 
iDiemand et ^20081 ISpringel et aljl2008btlKravtsovl2009b . 

However, these estimates of substructure abundance at dif- 
ferent locations within a halo, from the solar neighborhood to 
the outskirts of the Milky Way, are based on simulations that 
include only dark mat ter and do not acco unt for ordinary bary- 
onic material (but see lDolag et al.l l2008). This means that any 
dynamical coupling between dark matter substructures and 
the luminous components of galaxies has been ignored, even 
though this may alter the substructure abundance. A precise 
understanding of the impact of baryons on the substructure 
abundance is important to correctly assess the relevance of 
substructures for the evolution of galaxies, and for attempts 
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to detect dark matter through gravitational lensing or annihi- 
lation radiation. 

For example, recent studies suggest that cosmological in- 
fall of dark satellites through the disk of a typical galaxy 
could induce bars an d multi-armed spiral patterns in disks 
dDubinski et alj|2008l) . According to these simulations, the 
main agent producing transient features would be satellite 
passages through the inner part of the disk. Indeed, it is spec- 
ulated that the passage of a satellite through the disk induces 
a localized disturbance tha t could grow by the mechanism 
termed swing amplification (lToomrell981l) . Because the tidal 
effects of the satellites are generally small, this process is dis- 
tinct from interactions t hought to be responsib le for grand- 
design spirals like M51 ( Kaz antzidis et a l. 2009). The mech- 
anism proposed by Dub inski et aL[ (2008) i s closer in sp irit 
to the original ideas suggested by Julian & Toomre ( 1966) in 
which an object within a disk excites a response that is sub- 
sequently amplified. Unlike the interactions involving from 
large visible satellite galaxies, we are more interested in the 
process produced by the continuing bombardment of a galac- 
tic disk by the dark matter halo inhomogeneities that will keep 
on exciting ragged spiral structures. 

Here, we investigate the coupling between dark matter sub- 
structures and the luminous components of galaxies. We con- 
sider tidal processes that lead to a depletion of the substruc- 
tures orbiting in the inner regions of a galaxy by the presence 
of a disk and we discuss its astrophysical implications. 

Naively, one might expect that the luminous components of 
galaxies should play only a minor role in affecting substruc- 
tures because most of the mass is in the dark matter. How- 
ever, some gravitational interactions between bodies, espe- 
cially those involving tidal forces, depend not on the relative 
masses of objects, but on their relative densities. Whereas the 
dark matter is thought to be nearly collisionless, baryons can 
dissipate energy and angular momentum and collapse to high 
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densities. This must have been the case for the baryons in the 
disk of the Milky Way, which is far more concentrated in its 
inner parts than the dark matter that makes up the surrounding 
smooth halo. Thus, it is plausible that the baryons could have 
an impact on the dark matter substructures despite their small 
contribution to the overall matter density of the Universe, mo- 
tivating our investigation. 

Specifically, we carry out a number of fiducial numerical 
experiments in which we adiabatically grow a disk in a high- 
resolution collisionless simulation of a Milky-Way sized halo. 
This enables us to examine the impact of e.g. disk shocking 
on dark matter substructures that pass through the disk. We 
also follow the orbits of substructures in the high-resolution 
dark matter simulations of the Aquarius project, allowing us 
to identify the parameters and frequency of encounters with 
an assumed disk in the halo. This provides an independent 
check of our isolated disk-growth simulations, and also en- 
ables an investigation of the consequences of different mod- 
eling assumptions. 

In S|2] we discuss the numerical set-up of our disk-growth 
simulations, and in Sj3]we present the results. In Sj4]we com- 
pare with analytic estimates of the disk shocking effect based 
on tracking satellite orbits in cosmological simulation of halo 
formation. Finally, we discuss the implications of our results 
ing] 

2. NUMERICAL METHODS 

Our simulations were performed with GADGET3, a par- 
allel TreePM code developed to compute the evolution of 
dissipationless dark matter systems. A deta iled description 
of the code is avai lable in the literature (Springel 120051: 
Sprin gel et al.ll2008ab . Here, we note its essential features. 

GADGET3 is a cosmological code in which the gravitational 
field on large scales is calculated with a particle-mesh (PM) 
algorithm, while the short-range forces are provided by a tree- 
based hierarchical multipole expansion, such that an accurate 
and fast gravitational solver results. The scheme combines the 
high spatial resoluti on and relative insensitivity to clustering 
of tree algorithms dBouchet & HernquistHl988l) with the un- 
matched speed and accuracy of the PM method to calculate 
the long range gravitational field. Pairwise particle interac- 
tions are softened with a spline of scalelength h s , so that they 
are strictly Newtonian fo r particles separated by more than h s 
dHinquist & Katz 1989). For our applications, the gravita- 
tional softening length is kept fixed in comoving coordinates 
throughout the evolution of our halos. 

2.1. Setting the initial conditions 

Presently, galaxy formation is not well-understood and so 
it is not possible to perform an ab initio cosmological simu- 
lation that combines the evolution of dark matter and baryons 
and produces an object that closely resembles the Milky Way. 
For this reason, we instead employ simulations in which we 
follow the growth and subsequent gravitational collapse and 
virialization of dark matter halos self-consistently, but include 
the gravitational contribution from luminous galaxies in an 
idealized manner. In particular, we represent the baryons of a 
galaxy as a fixed disk potential, added to a dark matter simula- 
tion with parameters and temporal behavior chosen to match 
observational and theoretical constraints. This enables us to 
study the dynamical coupling between luminous galaxies and 
substructures in their dark matter halos. 

In our procedure, we first evolve cosmological simulations 
with dark matter only and locate halos with properties at the 



present day similar to those of the Milky Way. The halo de- 
scribed in detail below was originally selected within a low 
resolution, dark matter only simulation run in a concordance, 
A-dominated universe with parameters Q.„, = 0.27, A = 0.73, 
h = 0.7, Og = 0.8, and Clt, — 0.044. These cosmological 
parameters are consistent with the current set of constraints 
within their uncertainties, i n particular those from the WMAP 
1- and 5-year data analysis (Komatsu et al. 2009). 

The large dynamic range involved in cosmological simula- 
tions aimed at resolving the scales of halos like the one pre- 
sented here calls for techniques that concentrate the compu- 
tational power on the object of interest. This is achieved by 
a zoom-technique, where a large-scale simulation is done at 
low resolution and the relevant regions are identified within 
it, e.g. a Milky Way-sized halo such as in this work. The sim- 
ulation can then be repeated with much higher resolution in 
this region. 

The size of the box we employed, 10/z~ ! Mpc on a side, 
is quite small and starts to become non-linear at the present 
epoch, but it is still sufficient to provide realistic torques 
and accurate formation histories of Milky-Way sized halos. 
The virial mass of the halo we selected at z = to be re- 
simulated at higher resolution is 5.5 x 10 11 h~ l M & (the virial 
mass is here measured within a radius /? v ir ~ 160/i _1 kpc, 
enclosing an overdensit y of ~ 96 times the critical density 
dBrvan & Normanlll998l) ). 

Dark matter particle masses in the high resolution re- 
gions are 5.5 x 10 5 /z~'M Q , and the force resolution, i.e. the 
gravitational softening, is 200/z _1 pc. In total there are ~ 
10 6 particles within the virial radius. With our choices 
for particle number and softening, the smallest subhalos re- 
solved have typical masses of ~ 10 7 M fi . The chosen halo 
has a merger history and spin parameter reasonably repre- 
sentative of the global p opulation of halos in equilibrium 
(D'Onghia & Navarro 2007). The halo was selected with the 
only criterion that the redshift of the last major merger was 
jRil.5 and that there are no halos of similar or larger mass 
within a few virial radii (a major merger is defined here as 
having a 4:1 mass ratio). The halo spin parameter is 0.04 at 
z = 0, close to the averag e value » 0.035 for halos in large 
cosmological simulations dBett et al.ll2007l) . 

2.2. Putting a disk inside the galactic halo 

We began the cosmological simulation at high resolution in 
the region of interest at an early time and evolved it forward 
to the present day (and beyond), saving the state of the system 
at many intermediate times. We extracted the galaxy halo at 
the present time from the cosmological simulation, including 
all particles within 500 hr 1 kpc of the galactic center and then 
repeated the simulations, but adding the potential from a disk 
to the halo starting at a specific time. We repeated this pro- 
cedure by varying both the manner in which the galaxy was 
included as well as its properties. 

We have modeled the disk potential according to (Kuzmin 
1956): 



^/fi 2 + ( fl + v?+F) 2 ' 

where M is the disk mass, a is the scale-length of the disk and 
is taken to be 4.55/;~ 1 kpc, and b is the disk height, assumed 
to be 260pc. With these choices, this po tential matches that of 
the Galactic disk (Johnst on et alj[l99 5). We have tested dif- 
ferent prescriptions to place the disk inside the halo: we either 
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FIG. 1. — Dark matter distribution within 60 kpc of the center of the dark matter halo in two simulations, each at three different times. The columns are the 
initial state (left column; identical in both cases), after 3.7 Gyrs (middle column), and after 5.5 Gyrs (right column). Top row is a simulation with no disk, bottom 
row is one where a disk with mass 0. 1 times the dark halo was included. For the experiment with a disk, the disk was grown adiabatically over the initial 2 Gyrs 
to its final mass. 



centered on the halo center of mass, or on the particle where 
the halo potential is a minimum. In the former method we 
adopted an iterative method to locate the maximum density 
and centered the disk there. The different procedures gave 
similar results; for the final runs analyzed in this paper we 
adopted the minimum potential approach. 

In the simulations with disks, we grew the galaxy slowly 
(adiabatically) for 2 Gyrs to its full mass to avoid initial tran- 
sients. In our fiducial case, the disk grows to a mass of l/10th 
of that of the halo within its virial radius, but we also per- 
formed different simulations altering the properties of the lu- 
minous galaxy. In particular, we varied the mass of the disk 
relative to that of the dark matter halo, its radial extent, and 
the growth timescale of the disk. In general, we see the heat- 
ing processes becoming stronger as the mass of the disk or its 
scale length are increased, consistent with analytic expecta- 
tions (see below). 

We investigate the abundance of dark matter substruc- 
tures inside the halo as measured by the SUBFIND algorithm 
( Sprin gel et al.l200lT) . All our substructures consist of particle 
groups that are gravitationally self-bound and are overdense 
with respect to the local background. Every simulation parti- 
cle can be part of only one subhalo. We count substructures 
down to a minimum of 20 bound particles. 

3. RESULTS FOR DISK-GROWTH SIMULATIONS 

The generic outcome of our experiments is illustrated in 
Figure Q] which shows the appearance of halos at various 
times for cases in which we ignored the contribution from a 
disk (top row) and included it (bottom row). From left to right, 
we show the dark matter distribution initially and then after 
3.7 and 5.5Gyrs. The simulations shown in this figure as- 
sumed a disk with a final mass of 5.5 x lO 10 /!" 1 M n , i.e. 10% 
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FIG. 2. — Cumulative number of subhalos as a function of distance from 
the center of the host halo, with and without the disk. 

of the dark mass within the virial radius, appropriate for the 
Milky Way. 

As is visually apparent in FigureQ] the dark matter distribu- 
tion is significantly influenced by a disk. Within 2 Gyrs after 
reaching its final mass, the disk is able to erase most of the 
substructures inside the inner few tens of kpc of the halo. Be- 
cause this timescale is short compared to the age of the Uni- 
verse and the evolutionary timescales of real galactic disks, 
we expect that few, if any, low mass substructures would actu- 
ally be orbiting the innermost regions of a large spiral galaxy. 

This is made more explicit in Figure [2] where we compare 
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FIG. 3. — Cumulative number of subhalos above a given mass limit, for the 
evolved simulations with and without a disk within two virial radii. At late 
times, the disk can reduce the substructure abundance by nearly a factor of 
two. 




FIG. 4. — Spatial distribution of substructures within 150/i~'kpc of the 
galaxy center, comparing the case when a disk is not included (top set of pan- 
els) with the case where disk is included (bottom set of panels). For both, 
we show the distribution at the initial time t = (top rows) and after 3.7 Gyrs 
and 5.5 Gyrs, respectively (middle and bottom rows). Circle areas are pro- 
portional to the substructure masses. 

the radial cumulative abundance of substructures between the 
runs with and without disks. As time goes by, the depression 
of the subhalo function in the disk case increases, in propor- 
tion to the larger number of disk passages that have occurred. 
This substructure depletion is particularly strong in the inner 
parts of the halo, even though an effect is also noticeable in 
the outer parts. In Figure [3j we plot instead the cumulative 
substructure count as a function of mass, comparing again the 
simulations with and without a disk at two different times af- 
ter the start of the simulations. There is up to a factor of ~ 2 
reduction in the substructure abundance, and the effect is ap- 
proximately independent of substructure mass if we discard 
the measurement for the most massive bin, which is popu- 
lated by only one object and therefore allows no statistical 
conclusions. 

The depletion of the substructure abundance is also re- 
flected in the spatial distribution of substructures, which we 
show in Figure@]at different times for the cases without a disk 
(top) and when a disk with a final mass equal to 10% of the 
final halo mass included (bottom). The areas of the symbols 
are proportional to satellite mass. 

Figures|2][3]and|4]confirm that a disk accelerates mass loss 
by satellites, altering the subhalo mass function. Overall, the 



FIG. 5. — Energy of a satellite of 10 7 Mq as compared to its binding energy 
being affected by disk shocking produced by an exponential disk (blue dashed 
line), by a Kuzmin disk as in our simulations (solid red line), and by tidal 
shocking from the halo (green solid line) for different pericentric distances of 
the satellite orbit. 




1 10 100 

R [kpc/h] 



FIG. 6. — Radial dependence of the intermediate to major axis (b/a) (bot- 
tom panel) and minor to major axis (c/ a) (top panel) ratios of ellipsoids fitted 
to the halo when a disk is included. 

effect of the disk is to reduce the number of substructures and 
the masses of those that survive already after 3.7 Gyrs. More- 
over, with time the disk also compresses the dark matter dis- 
tribution in the center, further contributing to subhalo heating 
and accelerating the mass loss of these systems. A measure- 
ment of the halo radial dark matter density profile showed that 
the halo contracts and becomes denser by a factor of 2 within 
the inner regions owing to the gravitational potential of the 
disk. We also note that the presence of a disk makes the inner 
parts of the dark matter halo rounder, an effect that we further 
quantify below. 

We can identify three physical mechanisms which cause 
substructures to lose mass and to eventually be disrupted. 
First, as these objects orbit within the halo-galaxy sys- 
tem, they will be tidally stripped. The contribution 
of the halo to this process has been included in previ- 
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ous studies dHavashi et all 120031: iKazantzidis et all [2 004b; 
iRead et al.ll2006UPenarrubia et alJl2008t iMufloz et alJl2008f) 
but generally not tha t from the luminous galaxy (but see 
iTavlor & B abul 2001). Moreover, the adiabatic contraction 
of the halo owing to the disk enhances this process, an effect 
that is not accounted for in dark matter only simulations. 

The impact of this effect can be estimated crudely by as- 
suming that the orbiting substructure is truncated at some tidal 
radius r, where the differential tidal force of the halo-galaxy 
system is equal to the gravitational attraction of the satellite. 
For non-circular orbits the tidal radius is such that the av- 
erage density of the satellite within it is proportional to the 
mean density of the Galaxy at R p , the pericentric distance: 
Psub('"r) x Pgai(^p)- N-body simulations confi rm that subha- 
los are truncated at the radius given by r, (e.g. Sprin gelet al.l 
2008a), but we expect that this process will be dominant only 
for satellites orbiting with large pericentric distances. 

Second, the tidal field of the dark matter halo causes tidal 
heating of the particles in a satellite galaxy. Thi s process is 
most efficient for satellites on eccentric orbits (Gned inet alJ 
1999), which are the majority according to numerical sim- 
ulations, whereas it vanishes for systems on circular orbits. 
Within ~ 20kpc of the Milky Way center this process destroys 
substructures more rapidly than tidal stripping, owing to the 
strong tides that occur during close passages of the central 
dark matter cusp. We hence may also call this process halo 
shocking. 

Third, when an object plunges through the dense baryonic 
disk, it will be subject to a process known as disk shocking 
(Ostri ker et al.lll972l : lBinney & Tremaind 19871) which can, in 
some circumsta nces, produce a muc h stronger response than 
tidal stripping (Taylor & Babul 2001). This process is not ac- 
counted for in dark matter only simulations, and its relevance 
for the evolution of the substructure abundance has largely 
been ignored thus far. We will hence focus on this process in 
what follows. 

Whereas a slowly varying tidal field can strip loosely bound 
material through tidal stripping, a rapidly varying gravita- 
tional field, such as that arising when a body passes through 
a galactic disk, will induce gravitational shocks that add en- 
ergy to an object, changing its structure and accelerating mass 
loss. For an object orbiting in the inner regions of a galaxy on 
an orbit inclined with respect to the disk, disk shocking may 
indeed dominate the mass loss. This process i s known to in- 
fluenc e the structure of globular star clusters (Ostrik er et alJ 
fT972h . and can be more significant for dark matter substruc- 
tures. In the solar neighborhood, the galactic disk has a mass 
density p ~ 0.2 M Q pc~ 3 and the vertical scale-height with 
which the density decreases with distance above or bel ow the 
midplane is roughly 350pc (Bin ney & Tremainelll987l) . A 
clump of dark matter with mass of Msj 10 7 M Q and half- 
mass radius rwj ~ lkpc has a mean density within rw 2 of 

order p ~ 0.002 M Q pc -3 , 1000 times less than that of either 
the Milky Way disk locally or a globular cluster with M« 10 6 
M and r\n ~ lOpc. Therefore, disk shocking can be more 
important for the dynamics of dark matter substructures than 
for globular clusters. 

Indeed, this suggestion is supported by an additional con- 
sideration. The characteristic internal velociti es for a bound , 
virialized system are (v 2 ) ~ 0.4GM/r 1/2 (ISpitzeri 1 19691: 

Binney & Tremaine 1987), typically ~ 7.5kms~ 1 for a glob- 
ular cluster, and only 1.4kms _1 for a clump of dark mat- 
ter of this mass. When an object passes through a disk 



of vertical extent 2z, its constituents will move a distance 
Ar ~ 2zy/ (v 2 ) /Vj_ along their orbits, where Vj_ ~ 170kms _1 
is the mean vertical speed for typical inclined orbits t hrough 
the disk of the Milky Way dBinney & Tremain9 [l987). Stars 
in a globular cluster will thus move a distance Ar ~ 50 pc dur- 
ing a disk crossing, implying that tightl y bound stars will b e 
protected by adiabatic invariance (Binnev & Tremaine : 19871) . 
and that only the outer, more loosely bound material will be 
affected. However, because the internal velocities character- 
istic of mass clumps of dark matter are lower, the particles 
within these objects will move smaller distances along their 
orbits, Ar ~ 6 pc. This is much less than the half-mass radii 
of these substructures; therefore energy can be efficiently de- 
posited into the bulk of the mass of such an object by disk 
shocking, in some cases leading to its complete disruption. 

A substructure will acquire from a disk passage an energy 
per unit mass of order AE ~ 2(r 2 )g 2 /3V 2 , where g z is the 
local gravitational acceleration above and below the disk. We 
note that the second moment (r 2 ) of the radial particle dis- 
tribution of a subhalo is tightly correlated with its half-mass 
radius r\ji\ from the Aquarius N-body simulations we find 

(r 2 ) ~ (1.84r 1/2 ) 2 . The heating energy AE is to be com- 
pared with the total gravitational binding energy of the ob- 
ject, which we can for example estimate accordi ng to Ef, sa 
— 0.28GM 2 /ri/2i by fitting the density profile to a lHernquisl 
( 1990) model, although this choice makes little difference to 
our argument. An alternative is to assume virial equilibrium 
and estimate the binding energy from the total kinetic energy 
in random motions of the bound subhalo, i.e. Eh w M (v 2 ) /2. 

For a substructure of mass of 1O 7 M , we find AE ~ E b 
per orbit at 10 kpc from the center. If the disk mass den- 
sity declines exponentially with radius, then a corresponding 
estimate gives AE w Q.lEj, per orbit at a distance of 20 kpc 
from the center. Thus, complete disruption from disk shock- 
ing alone could occur at or within the solar radius during a sin- 
gle orbit, whereas w 10 orbits would be required to disrupt a 
substructure orbiting at a distance of 20 kpc. For objects pass- 
ing through the disk in the solar vicinity on modestly inclined 
or eccentric orbits, the azimuthal period is TJj, ~ 2 X 10 8 years, 
implying a disruption timescale of order ~2x 10 9 years. This 
should be regarded as an upper limit because tidal stripping 
will accelerate mass loss when the substructure is above or be- 
low the disk, and because this estimate ignores the reduction 
in binding energy following each disk shock. These calcula- 
tions may overestimate the rate of depletion for the subhalos 
with the half mass radius larger than the disk thickness. In 
this case the maximum amount of energy a substructure can 
receive per unit mass will be limited by the region of the sub- 
structure that is compressed/shocked at a given time which is 
of the order of disk thickness. Thus, the overall effect of disk 
shocking for subhalos larger than 10 7 might be overestimated 
in this model. 

A similar estimate applied to the more massive dark matter 
substructures in our simulations suggests that the disruption 
timescale for these objects is also short near the disk, explain- 
ing the rapid depletion of substructures seen in the inner re- 
gions of Figure[U Indeed, in the simulation shown in Figure[T] 
with a disk, only one subhalo with maximum circular veloc- 
ity larger than 10km s -1 survives within 30 kpc of the center 
of the galaxy. As it happens, this is in agreement with the 
one known Milky Way satellite on such a tightly bound orbit, 
Sagittarius at 24 ± 2 kpc. Thus, it does not appear that there is 
an inner missing satellite problem in the Milky Way. 
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FIG. 7. — Examples of radial orbits for four surviving substructures in the Aq-A-2 simulation of Springel et al. i2008a). In each panel, the solid black line shows 
the physical distance of the substructure center to the center of the main progenitor halo of the forming Milky- Way sized halo. The virial radius of this growing 
main halo is shown as a blue line. The red circles mark the output times used in the analysis; the actual orbit is constructed with an interpolation technique based 
on these output times. The thin vertical lines mark the pericentric passages that are identified by our orbit tracking technique. 

In Figure|5]we show the efficiency of the different processes 
causing mass loss and destruction of satellites as a function of 
the pericentric distances of their orbits. For a substructure or- 
biting within 20kpc of the galactic center the halo shocking 
process is stronger than tidal stripping. However, at slightly 
smaller radius, disk shocking becomes the dominant destruc- 
tion process overall. Because of the presence of a disk in 
the inner regions, the Galactic halo contracts and becomes 
slightly denser in the inner regions, further enhancing the im- 
pact of tidal heating and stripping, although we note that our 
estimates in Figure[5]do not account for this. 

3.1. Halo shape 

Finally, we also note that the dark halo in Figure Q] ap- 
pears to be more spherical when the disk is included (in par- 
ticular, the middle and third columns of the bottom row) as 
compared to the simulation without a disk (top row), in qual- 
itative agreement with previous studies which have tried to 
include baryonic material in live halos (Springel et al.l 120041: 



iKazantzidis et al . 2004a; Maccio' et al.ll2007t iDebattista et all 
2008; Aba di et all 120091) . Indeed, the halo responds to the 
presence of the disk by becoming significantly more spheri- 
cal. This is quantified in Figure|6]where we measure the radial 
dependence of axial ratios of ellipsoidal surfaces of the halo 
as a function from the halo center when a disk is included. 
The panels in Figure [6] display the intermediate to major axis 
(b/a) (bottom panel) and minor to major axis ratio (c/a) (top 
panel) of the halo measured at different distances from the 
center. The presence of the disk turns the halo from a triax- 
ial to more spherical object within 30 — 40kpc of the center, 
already after 3.7Gyrs (the disk is fully in place after 2 Gyrs). 
The halo shape is not independent of radius. Indeed, this halo 
becomes more round with a baryonic component only within 
30 — 40kpc, roughly the size of the disk and it is not affected 
as far out the virial radius. Hence, our model suggests that 
halos around disk galaxies should be round only within the 
regions occupied by the disk and not necessarily far beyond 
that. 
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4. DISK PASSAGES IN COSMOLOGICAL SIMULATIONS 

In order to extend the dynamic range of the masses of sub- 
structures for which we can study the disk-shocking effect, 
we follow the orbits of subhalos in one of the high-resolution 
numerical realizations o f Milky Way sized halos calculated 
in the 'Aquarius' project (Springel et al. 2008b HNavarro et aT] 
2008). In this suite of dark matter only simulations, extensive 
resolution studies for 6 different halos of mass ~ 10 12 M fi 
were carried out, dark matter halos and substructures were 
found at a large number of output times, and detailed merger 
history trees for all these halos and subhalos were constructed. 
Using these trees allows us to accurately follow the orbits of 
individual subhalos and to determine how often they are ex- 
pected to cross a fiducial disk assumed to be present in the 
galaxy-sized halos. 

We here focus most of our analysis on the 'Aq-A-2' simula- 
tion, which has a mass resolution of 10 4 h~ l M&, and a virial 
mass ofM 20 o = 1.34 x IO^/t'Mq within R 2 oo = 180/^'kpc 
(the radius enclosing a mean overdensity of 200 relative to the 
critical density) for its final Milky Way-sized halo. Based on 
the substructure catalogs and merger trees, we construct con- 
tinuous orbits for all substructures by interpolating the sub- 
halo trajectories with the discrete set of coordinates and ve- 
locities available for the output times. We use reconstructed 
orbits based on 128 outputs between redshifts z ~ 20 and z — 
for the analysis shown here, but we have checked that using 
up to 1024 outputs instead does not change the estimated disk 
encounter parameters in any significant way; i.e. the time res- 
olution of 128 outputs is sufficient to obtain sufficiently accu- 
rate orbital tracks for the purposes of this analysis. 

In Figure|7] we show examples of radial orbits for four typi- 
cal substructures with mass of around 10 8 h~ 1 M & that end up 
within the virial radius of the final halo. We clearly see that 
the subhalos make several orbits, sometimes reaching quite 
small pericentric distances to the halo center. Note that even 
substructures that presently have a comparatively large dis- 
tance from the center may have had a close encounter with 
the halo center at some earlier time. For example, the sub- 
halo shown in the bottom right panel has now a distance of 
r ~ 90/z~ ! kpc, but it passed by the center at a distance of 
less than 2/z~ 1 kpc as early as z — 2.5. Such close passages 
at small distances to the center may then lead not only to sig- 
nificant halo shocking (which the substructures shown in this 
plot in fact did survive) but also to disk shocking. 

In order to obtain an estimate for the cumulative effect of 
the latter, we examine the orbits of all substructures that end 
up inside the virial radius at z = and which have a mass 
above M su b = 10 6 /!~ 1 M Q at the final time. For each passage 
through a fiducial plane containing the disk of the main halo 
(say the xy-plane, but choosing a different orientation gives 
consistent results), we register the distance R to the center, the 
velocity V± perpendicular to the disk, the subhalo's current 
half-mass radius r^n and its squared 3D velocity dispersion 

(v 2 ). This allows us to characterize the strength of the disk- 
shocking as 

AE _ (1.84r 1/2 ) 2 [47tG£(fl)] 2 

E 3<v 2 }V 2 ' { ' 

where we assumed the infinite sheet approximation to relate 
the disk's surface mass density T.(R) to the vertical gravita- 
tional field above and below the disk. For systems that expe- 
rience multiple disk transitions, we simply add up the differ- 
ent AE/E values. A significant disk shocking effect can be 



expected when AE/E reaches values of order unity. We shall 
assume that AE/E > 1 implies certain disruption, but we note 
that also smaller values, AE/E > 0.1, should lead to signifi- 
cantly accelerated mass loss and earlier destruction. 

The distribution of AE/E values obtained in this way is 
broad, with many subhalos having very small values, sug- 
gesting that for them disk shocking is unimportant. However, 
we find that about 14.1% of the subhalos with mass above 
M sub = 10 6 ^ 1 M Q have AE/E larger than 1.0, and 24.2% 
have a value larger than 0. 1 . 

In Figure[8] we show the estimated impact of disk-shocking 
in terms of the cumulative radial count of subhalos, compar- 
ing the counts of all subhalos to those where systems with 
AE/E> l.OoiAE/E > 0.1 are omitted. Interestingly, the cu- 
mulative effect is clearly stronger in the inner than in the outer 
parts. In fact, for the small mass subhalos that dominate the 
M su b > h 1 1O 6 M sample, the reduction at radii of ~ 30kpc 
reaches a factor of ~ 3 — 4, whereas it is only about 25% at 
the virial radius. We also include results restricted to the more 
massive subhalos M su b = lO 8 /z _1 M in Figure[8] This shows 
a qualitatively similar trend, but with poorer statistics. We find 
that the more massive subhalos tend to be less affected by disk 
shocking, which can be understood based on their more recent 
infall time and hence smaller number of passages through the 
disk. 

We note that the precise strength of the disk shocking ef- 
fect depends, of course, on the adopted disk model. In the 
results above, we have assumed a fixed exponential disk with 
scale radius 5 hr 'kpc and with mass equal to 10% of the virial 
mass of the final halo. Alternatively, we have explored mod- 
els where we assumed that the disk mass grows in proportion 
to the virial mass of the halo, and the disk size scales with the 
halo's virial radius. This moderately reduces the strength of 
the effect; we then find that 17.1% instead of 24.2% of the 
subhalos reach a value AE/E > 0.1. 

Another very interesting question is whether disk shocking 
in progenitor systems other than those of the main halo may 
be important for the final substructure population. Note that 
the results above have only estimated disk shocking from the 
central disk, but perhaps many substructures fall in as part of 
groups, in which case they could already have experienced 
severe tidal shocks owing to the baryonic structures in these 
groups. 

To investigate this idea, we first ask the question of what 
is the fraction of substructures that have fallen in as a sub- 
halo within a larger system, as opposed to being accreted as 
a main halo on their own. This can be readily addressed with 
the merger trees of the Aquarius halos. In Figure|9] we show 
the fraction of subhalos that were accreted as a substructure of 
an infalling group, as a function of their present day mass for 
all subhalos in the virial radius. We compare results for sim- 
ulations of the same halo at different numerical resolutions, 
finding good convergence. We see that the fraction of sub- 
structures that come in already as a subhalo in a larger group 
is actually quite small, only about 10-20% over the mass range 
considered, with a slight trend to i ncrease towards smaller 
masses (see also lLudlow et al.ll2009l) . In hindsight this is per- 
haps not too surprising, because for any given mass, there are 
always more main halos as genuine subhalos, and the merger 
hierarchy i n cold dark matter mo dels is actually surprisingly 
'shallow' dAngulo & Whitdl2009t) . 

However, this also means that disk shocking in progeni- 
tor systems different from the progenitors of the main halo 
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(which form the 'main stem' of the merger tree) is unlikely to 
be very important. 

We explicitly confirm this by extending our estimates of 
disk shocking to include all progenitor groups and not only 
the main halo. To this end, we extract disk passage parame- 
ters also when the reconstructed orbit of a subhalo passes by 
the center of a main halo different from the main progenitor 
of the primary galaxy. For simplicity, we assume that there 
are also disks present in these secondary progenitor systems, 
and that they always have a disk mass equal to 10% of their 
current virial mass, with a size scaled in proportion to their 
virial radius. We then find that only a relatively small number 
of additional subhalos are predicted to be strongly affected by 
disk shocking. In particular, when using the AE/E > 0.1 cri- 
terion, the number of affected subhalos goes up from 24.2% 
to 24.8% when the disk in the main halo is kept fixed, while 
it goes up from 17.1% to 17.8% when the disk is assumed to 
always scale with halo size. 

Comparing the results from this analysis of substructure or- 
bits in the Aquarius halos to our direct disk-growth simula- 
tions, we find in general reasonable agreement. In particular, 
the effect estimated from the Aquarius simulations for the in- 
ner regions of the halo is of comparable magnitude to that 
found in our disk-growth experiments, even though it appears 
a bit weaker overall. But this is not unexpected. For one, the 
disk-growth simulations experience a central compression of 
the dark matter halo which causes an additional increase in the 
efficiency of halo shocking, an effect that we have not taken 
into account in the analysis in the present section. In addition, 
disk shocking will accelerate the mass loss of surviving sub- 
halos, which is also an effect that we have not included here. 
We emphasize that the calculations in this section are based 
on analytic estimates and involve several approximations. A 
self-consistent analysis of the orbital properties of subhalos 
and their mass loss using numerical simulations when the disk 
is included in a Milky Way halo is required to study these is- 
sues in more detail. 

5. DISCUSSION AND CONCLUSION 

Tidal halo shocking and disk shocking are dominant effects 
compared to tidal stripping for fairly eccentric orbits and have 
not been included in all previous simulations of dark satellites 
orbiting luminous galaxies like the Milky Way. Our work 
shows that these processes can be efficient at depleting the 
satellite population within ~30 kpc of the Milky Way center. 

While the effects we model in our simulations can extend 
further than the radial scale-lengths of disks, because many 
substructures are on moderately radial orbits, our calculations 
do not resolve the "missing satellite problem" on larger spatial 
scales. In principle, the processes we have investigated here 
could propagate from smaller to larger mass scales if one ac- 
counted for visible galaxies in the substructures that merged to 
form the Mily Way halo, driving pre-destruction of substruc- 
tures within subhalos and leading to a stronger depletion of 
satellites in the outer regions of large halos. However, our es- 
timates indicate that most substructures are accreted as main 
halos and not as subhalos, limiting the overall importance of 
this effect. We do note, however, that if luminous satellite 
galaxies in the Milky Way formed more efficiently in sub- 
ha loes accreted as pa r t of a Magellanic Group, as proposed 
by D'Onghia & Lakd (120081) . then it would not be surprising 
that most substructures would remain dark, even if they are 
not destroyed. 

Our findings have interesting implications for understand- 
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FIG. 8. — Cumulative radial distribution of satellites with mass larger than 
10 6 /i _1 M Q (black lines) and larger than 10 8 /r 'M Q (red lines) in the Aq-A-2 
high-resolution simulation of the Aquarius project. The cumulative count is 
shown for the original subhalo population (solid lines), and when subhalos 
are excluded for which the estimated disk shocking on their orbits exceeds 
AE/E = 1.0 (dashed) and AE/E = 0.1 (dot-dashed), respectively. 



0.1 



I 1 1 


1 1 1 


i 1 _ 


; Aq-A-2 






- Aq-A-3 






- Aq-A-4 

1 i i 


i i i 


i i 



I 

10 5 10 6 10 7 10 8 10 9 

M iub [ M e ] 



FIG. 9. — Fraction of substructures present in the main halo at z = that 
were accreted as subhalos of another infalling structure, as opposed to being 
accreted as a main halo. We compare results at different numerical resolu- 
tions. For runs Aq-A-3 and Aq-A-4, the mass resolutions are factors of 3.6 
and 30 times worse than in our default simulation Aq-A-2, respectively. The 
dashed-line marks a power-law fit to the results for Aq-A-2. 

ing the nature of spiral structure in disk galaxies. Recent 
simulations have shown that dark satellites impacting a disk 
could play an important role in generating the multi-armed, 
global spirals seen in disk galaxies. Our work demonstrates 
that the majority of the satellites are depleted in a short time 
when they pass through the disk. This raises the question of 
whether an early impact of satellites inside, or later impacts 
outside, the disk could still provide a source of perturbations 
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capable of maintaining global spiral patterns in disk galaxies, 
or whether a different source is required. Interestingly, an al- 
ternative source might be molecular gas clumps with masses 
of 10 6 - 1O 7 M which are preferentially found close to the 
spiral arms in the Mi lky Way or infalling ha lo clouds which 
have similar masses dKeres & Hernquislll2009l) . 

Although we do not account for the presence of live disks 
we do not expect the passages of the tightly bound substruc- 
tures through the disk before these objects ar e destroyed to 
have a significant effect on disk kinematics. Hopki ns et al.1 
(2008) showed that the heating done to disks by smaller per- 
turbers is second order in the mass ratio, so heating of the 
disk by substructures in the mass range we have explored is 
likely to be unimportant, as also shown by Velazqu ez & Whitel 
( 1999) using N-body simulations. 

Another implication of our work is that the depletion of 
clumps within a few dozen of kiloparsecs from the Milky Way 
center might lead to a revision of the mass fraction embedded 
in substructures as implied by gravitational lensing. These ob- 
servations have been done for early-type galaxies where there 
is no disk and might lead to different values when applied 
to late-type/disk galaxies. Furthermore, recent work estimat- 
ing the flux ratio anomalies in gravitational lensing owing to 
the presence of substructure points out that the inner subhalo 
abundance predicted by CDM-only simulations app ears actu- 
ally s lightly too low to account for these anomalies dXu et alj 
2009). As we have shown here, because baryonic processes 
can reduce the inner substructure abundance significantly, it 
might turn out that this problem is even more severe. On the 
other hand, there are cases where these anomalous flux ratios 
are likely the result of visible substructures. In particular, in 
CLASS B2045 +265, one of the mo st extreme anomalous flux 
ratio systems. [McKean et al.l d2007) find a small dwarf galaxy 
companion which is sufficient to fully explain the observed 
flux ratios without the need for any additional dark substruc- 
tures; 

iDiemand et all d2008l) and ISpringel et ail (l2008bl) have 
computed the dark matter annihilation signal from the galac- 
tic substructures and from the smooth dark matter halo of 



the Milky Way using very high resolution cosmological Via 
Lactea II a nd Aquarius simulati ons, respectively. It was 
pointed out (Springel et al. 2008b) that the signal-to-noise for 
detection of dark matter annihilation from the galactic center 
is in general much larger than for any of the substructures, 
but in case substructures are detectable, their typical distance 
from the Solar circle is predicted to be around 20kpc. Hence 
this detectable substructure population would be strongly af- 
fected by the disc shocking effect. Other models of gamma- 
ray emission built on the Via Lactea II simulation predict that 
the number of dark clumps detectable from the Fermi-LAT 
is between 1 and 10 in five years of operation , for an opti- 
mistic dark matter scenario (iPieri et alJl2009bl) . but see also 
iKuhlen et all d2009l) for an even larger estimate assuming an 
additional strong boost of the annihilation rate due to unre- 
solved dark clumps. Our results reduce these values by a fac- 
tor of 3 for subhalos of mass 10 M , and by a factor of 2 
for larger masses. This strengtens the expectation that de- 
tection of the dark matter annihilation signal from the galac- 
tic center should be easier than from nearby substructures 
(Spri ngel et aTll2008bl: IPieri et al.ll2009al) . 

Ultimately disk shocking predicts that luminous satellite 
galaxies orbiting through galactic disks should leave long tails 
of stars and can be disrupted in a few Gyrs, a timescale much 
shorter than the tails live. This possibility is supported by re- 
cent observations of loops of tidal tails around NGC 5907, 
which may be t he remains of a ghost dwarf galaxy that has 
not been found dMartfnez-Delgado etai1l2008l) 
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